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ABSTRACT 

^ We present multiwavelength investigation of morphology, physical-environment, 

stellar contents and star formation activity in the vicinity of star-forming region Sh 
^ 2-100. It is found that the Sh 2-100 region contains seven H n regions of ultracompact 

and compact nature. The present estimation of distance for three H n regions, along 
^ with the kinematic distance for others, suggests that all of them belong to the same 

P-i molecular cloud complex. Using near-infrared photometry, we identified the most 

O probable ionizing sources of six H ii regions. Their approximate photometric spectral 

type estimates suggest that they are massive early-B to mid-0 zero-age-main-sequence 
, ^, stars and agree well with radio continuum observations at 1280 MHz, for sources 

^ whose emissions are optically thin at this frequency. The morphology of the complex 

> shows a non-uniform distribution of warm and hot dust, well mixed with the ionized 

^ gas, which correlates well with the variation of average visual extinction (~ 4.2 - 97 

mag) across the region. We estimated the physical parameters of ionized gas with the 
help of radio continuum observations. We detected an optically visible compact nebula 
2 located to the south of the 850 /im emission associated with one of the H u regions and 

Q the diagnostic of the optical emission line ratios gives electron density and electron 

!• temperature of ~ 0.67 x 10^ cm^^ and ~ 10"^ K, respectively. The physical parameters 

• '-j suggest that all the H n regions are m different stages of evolution, which correlate well 

% 
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with the probable ages in the range ~ 0.01 - 2 Myr of the ionizing sources. The spatial 
distribution of infrared excess stars, selected from near-infrared and IRAC color-color 
diagrams, correlates well with the association of gas and dust. The positions of infrared 
excess stars, ultracompact and compact H ii regions at the periphery of an H i shell, 
possibly created by a WR star, indicate that star formation in Sh 2-100 region might 
have been induced by an expanding H i shell. 

Subject headings: dust, extinction-galaxies: star clusters-H n regions-infrared: ISM- 
ISM: individual (SlOO)-radio continuum: ISM-stars: formation 



1. Introduction 

Massive star formation is poorly understood as compared to low-mass stars, because their 
formation takes place in the dense core of a molecular cloud of high visual extinction, usually 
observable at far-infrared (FIR) to millimeter wavelengths and their evolutionary time scales are 
much shorter (< 10^ yr) (e.g., Bemasconi & Maeder 1996). The radio free-free emission in terms 
of ultracompact H ii (UCH n ) regions, represents a later evolutionary sequence of a massive pro- 
tostar outlined by Beuther et al. (2007), where the associated high mass protostar may still be 
in the accretion phase (e.g. Shepherd & Churchwell 1996) or has already ceased the accretion 
(Garay & Lizano 1999). The UCH ii region can further evolve into the less obscure stage of com- 
pact H u (CH n ) and classical H n regions (Garay & Lizano 1999) by the disruption of associ- 
ated gas and dust, revealing both, the embedded high-mass and lower mass stellar population at 
shorter wavelengths. Most of the high-mass star-forming regions in the Galaxy lie at a distance of 
more than 1 kpc. As a result, the problem in interpreting observations of massive star formation 
originates in source confusion at the core of distant molecular clouds due to their cluster mode 
formation. Statistically, it has been shown that the most luminous protostars form in molecular 
clouds associated with evolved H ii regions (Dobashi et al. 2001), where the interplay between 
the stars' radiation field and associated gas and dust makes the environment even more complex. 
Therefore, a detailed study of the star-forming region hosting young massive stars, using optical, 
near-infrared (NIR), mid-infrared (MIR) and radio bands, is necessary to understand the high-mass 
star formation scenario in these complexes. 

In this paper, we have studied a young star- forming region (SFR) K3-50 shown in Fig. 1, 
which contains a group of H n regions, namely A, B, C, D, E, and F (cf. Israel 1976). The 
H n region C itself consists of two UCH n regions, CI and C2 (Harris 1975). The kinematic dis- 
tance towards K3-50 ranges from 7.9 to 9.3 kpc (Harris 1975; Bronfman et al. 1996; Araya et 
al. 2002). In the present study we have adopted a widely used value of distance of ~ 8.7 kpc 
for our analysis. We also derived distances of three H n regions that come out to be close to the 
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value of 8.7 kpc. The group of H ii regions is located within an area of radius ~ 3'. 5 centered on 

Qr2ooo = 20''01'"42^ ^2000 = +33°33'49", in the proximity of an evolved diffuse nebula Sh 2-100 
of ~ 4' in size (Sharpless 1959). Sh 2-100 is itself a part of the large molecular cloud complex 
W58, which consists of several classes of H n regions with widely varying physical parameters. 
The observed radio luminosity of the complex can be explained by the presence of OB association 
(cf. Israel 1976). In Fig. 1 one can notice that K3-50D, K3-50E, and K3-50F are associated with 
nebulosity. A faint nebulosity is also seen ~ 15.5" away in the southern direction of K3-50C2, 
whereas K3-50A appears as a stellar point-like source and is marked with a circle. The regions 
K3-50B and K3-50C1 are optically invisible radio sources, implying that the extinction is high 
enough to obscure the regions. By comparing the expected infrared Bm line fluxes as predicted 
from radio continuum fluxes, to the observed Bm line fluxes, Roelfsema et al. (1988) estimated 
visual extinction of 15, 26, 97 and 32 mag towards A (over 25" aperture), B (over 5" aperture), CI 
and C2 (over 15" aperture) respectively, while the extinction towards D was found to be 2 mag (5" 
aperture). Observations at radio wavelengths suggest that these regions harbor at least one massive 
OB star (Harris 1975; Israel 1976; DePree et al. 1994). The UCH n regions, astronomical masers, 
outflows, detection of high density tracer molecules and infrared (ER) sources with IR-excess are 
the trademarks of young star-forming regions. The evidence of weak bipolar molecular outflow by 
Phillips & Mampaso (1991) with CO / = 2 - 1 transition, ionized outflow by DePree et al. (1994) 
with H76 alpha radio recombination line, detection of molecular core mapped in the HCO"^ and 
CS (2-1) lines by Vogel & Welch (1983) and Bronfman et al. (1996) suggest the youthfulness of 
the Sh 2-100 region. These observations along with the presence of astronomical masers (H2O by 
Kurtz & Hofner 2005; OH by Baudry & Desmurs 2002) imply active star formation is going on 
around Sh 2-100. The current star formation activity in the vicinity of Sh 2-100 region is probably 
the eff^ect of feedback from the earlier generation stars of the W58 complex (Israel 1980). 

Though there are several studies in the radio and infrared domain, most of them are concen- 
trated on K3-50A and K3-50C. To continue our multiwavelength investigations of massive SFRs 
(cf. Ojha et al. 2004a; Ojha et al. 2004b; Ojha et al. 2004c; Tej et al. 2006; Samal et al. 
2007; Pandey et al. 2008), we have studied an area of 8' x 8' of Sh 2-100 region to understand 
the star formation activity in this region. This paper presents new results from optical and in- 
frared photometry, optical spectroscopy and low frequency radio continuum observations. Based 
on these observations, we have carried out a detailed study of the ionizing sources of individual 
H n regions, physical characteristics and the nature of associated lower mass population. With the 
following layout of the paper we tried to interpret the possible star formation scenario of W58 
cloud complex. In Sect. 2, we describe our observations and the reduction procedures. In Sect. 
3, we discuss other available datasets used in the present study. Sect. 4 describes the general 
morphology of the region. Sect. 5 describes the infrared and radio components associated with 
the region. In Sect. 6, we discuss the nature of individual regions. Sect. 7 is devoted to general 
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discussion and star formation scenario in W58 complex. We present the main conclusions of our 
results in Sect. 8. 



2. Observations and data reduction 

2.1. Optical photometry 

UBV CCD photometric observations were performed for the Sh 2-100 region on 2006 October 
24 and 25, using the 2K x 2K CCD system at the f/13 Cassegrain focus of the 104-cm Sampur- 
nanand telescope (ST) of ARIES, Nainital (India). The 0.37 arcsec pixel"' plate scale gives a field 
of view (FOV) of ~ 12'.6 x 12'. 6 on the sky. The observing conditions were photometric and the 
average FWHM during the observing period was ~ 1".7- 2".0. To improve the signal-to-noise 
(S/N) ratio, observations were made in 2 x 2 pixel binning mode. Several bias frames and twilight 
flat field exposures were taken during the observations. We observed the standard area SAlOl 
(Landolt 1992) on 2006 October 25 several times during the night to determine the atmospheric 
extinction coeflicients and to calibrate the CCD systems. The initial processing of the CCD im- 
ages was done using the IRAlj^ data reduction package. Then, for a given filter, frames of same 
exposure time were co-added to improve the S/N. Photometric measurements of the stars were 
performed using DAOPHOT II package of MIDASj^ The stellar images were well sampled and 
a variable point spread function (PSF) was applied from several uncontaminated stars present in 
the frames. The calibration from instrumental to standard system was done using the procedure 
outlined by Stetson (1987). The standardization residuals between the standard and transformed V 
magnitude and colors are less than 0.02 mag. We used only those stars having photometric error 
in magnitudes < 0.1 in the present analyses. 



2.2. Optical spectroscopy 

Using the Himalaya Faint Object Spectrograph Camera (HFOSC) available with the 2.01 m 
Himalayan Chandra Telescope (HCT), operated by the Indian Institute of Astrophysics, Bangalore 
(India), we obtained optical spectra of possible ionizing sources of K3-50D, K3-50B and optical 
nebula associated with K3-50A on 2007 June 18 and 2007 October 29. The instrument is equipped 
with a SITe 2K x 4K pixel CCD. The central 2K x 2K region is used for imaging which corre- 



'IRAF is distributed by the National Optical Astronomy Observatories, USA 
^MIDAS is developed and maintained by the European Southern Observatory 
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sponds to a FOV of ~ 10' x 10', with a pixel size of 0".3. The spectra were obtained with a grism 

o o , 

(Gr 7) in the wavelength range of 3500-7000 A at a spectral dispersion of 1.45 A pixel" . The ex- 
posure time was 900s for each source. The spectroscopic observations were carried out under good 
photometric sky conditions. The one-dimensional spectra were extracted from the bias-subtracted 
and flat- field corrected images using the optimal extraction method (Horne 1986) of IRAF. Wave- 
length calibration of the spectra was done using FeAr lamp source to identify emission lines of the 
lamp and we fitted a lower order polynomial to find a dispersion solution. For flux calibration of 
the nebular spectrum, we used spectrophotometric standard star BD 284211 (Massey et al. 1988) 
observed on the same night. 



2.3. Near-Infrared Observations 

Imaging observations oftheSh 2-100 region centered on tt2oo() = 20''01™44^ ^2000 = +33°33'07", 
in J (A= 1.25 yum), H (A = 1.63 fim), and = 2.14 /um) bands were obtained on 2008 July 21 
with the 1.4 m Infrared Survey Facility (IRSF) telescope and SIRIUS, a three-color simultaneous 
camera equipped with three 1024 x 1024 HgCdTe arrays. The FOV in each band is ~ 7'.8 x 7'.8, 
with a pixel scale of 0".45 at the f/10 Cassegrain focus. Further details of the instrument are given 
in Nagashima et al. (1999) and Nagayama et al. (2003). A sky field was also observed centered 
on orjooo = 20''00'"44^ ^2000 = +33°31'36" and used for sky subtraction. We obtained 120 dithered 
frames each of 10s, thus giving a total integration time of 1200s in each band. The observing con- 
ditions were photometric and the average FWHM during the observing period was ~ l".8-2".0. 
Data reduction was done using pipeline software based on IRAF package tasks. Dome flat-fielding 
and sky subtraction with a median sky frame were applied. Photometry of point sources was per- 
formed using the PSF algorithm ALLSTAR in the DAOPHOT package (Stetson 1987) within the 
IRAF environment. The PSF was determined from 15 relatively bright and isolated stars of the 
field. We used an aperture radius of 1 FWHM, with appropriate aperture corrections per band for 
the final photometry. 

For photometric caUbration, we used 54 isolated 2MASS Point Source Catalog (PSC) sources 
(see §3.1) with mag range 12-14 and having error < 0.05 mag. The same 2MASS sources were 
used for absolute position calibration and a position accuracy of better than +0".l was achieved. 
The sources with mag < 1 1 were saturated in our image. For such bright sources, 2MASS PSC 
data were used. For one source (i.e., a possible ionizing source of K3-50A), we have used the data 
(star id 68) from Alvarez et al. (2004), since the source was found to be « 1.3 mag and 2.2 mag 
brighter than /T-band and //-band, respectively, than in Alvarez et al. (2004). Alvarez et al. (2004) 
observation covers an area of ~ 80" x 80" in the vicinity of K3-50A with a resolution (FWHM) 
of 0".22. A comparison of present photometry with that of 2MASS photometry as a function of / 
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mag for the common sources is shown in Fig. 2, where the 2MASS sources of photometric error 
< 0.1 mag and 'read-flag' values between 1 to 3 are selected to ensure good quality. We find that 
the photometric scatter in magnitude increases at the fainter end, whereas the mean photometric 
uncertainties are < 0.1 in all the three bands. To plot the data in color-color (CC) and color- 
magnitude (CM) diagrams, the JHK^ data were transformed to California Institute of Technology 
(CIT) system using the relation given at the 2MASS website]^ 

2.4. Optical and NIR narrow-band imaging 

Optical CCD narrow-band images of the nebula around Sh 2-100 region were obtained in Ha 
{X = 6563 A, ^A = 100 A) and [S II] {A= 6724 A, ^A = 100 A) on 2007 June 18 with HFOSC on 
the HCT. We used i?-band as a nearby continuum for the continuum subtraction. 

Narrow-band NIR observations were carried out in Bry (A = 2.16 yum, AA = 0.019 //m) and K 
continuum (A = 2.264 /urn, AA = 0.054 /zm) with HCT on 2004 October 10, using the NIR imager 
(NIRCAM), which is a 0.8 - 2.5 yum camera with 512 x 512 HgCdTe array. For our observations, 
the NIRCAM was used in the Camera-B mode which has a FOV of ~ 3'. 6 x 3'. 6. We obtained 
several dithered (by 20") frames of the target centered on azooo = 20''01'"47', ^2000 = +33°31'38" 
in order to remove bad pixels, cosmic rays and to eliminate the presence of objects with extended 
emission for construction of sky images. Typical integration times per frame were 90s, and 20s 
in the Bry, and K continuum bands, respectively. The images were co-added to obtain the final 
image in each band. We also obtained several dithered sky frames close to the target position in 
both the bands. All images were bias and dark- subtracted and flat-field corrected. The continuum 
subtraction of the images was done after aligning and PSF matching. 

Details of the log of observations are given in Table 1. 



2.5. Radio continuum observations 

In order to probe the ionized gas component, radio continuum interferometric observations 
at 1280 MHz and 610 MHz were carried out on 2004 January 25 and 2004 April 07, using Giant 
Metrewave Radio Telescope (GMRT) array. The GMRT has a "Y"-shaped hybrid configuration of 
30 antennae, each 45 m in diameter. There are six antennae along each of the three arms (with arm 



^www.ipac.cal tech.edu/2mass/releases/allsky /doc 
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length of ~ 14 km). The rest of the twelve antennae are located in a random and compact 1x1 km^ 
arrangement near the center. Details of the GMRT antennae and their configurations can be found 
in Swarup et al. (1991). The largest angular scales to which the GMRT is sensitive are 8' and 17' 
at 1280 and 610 MHz, respectively. 

The flux and phase calibrators were observed in order to derive the phase and amplitude gains 
of the antennae and the details of the observations are given in Table. 2. 3C48 and 3C286 were 
used as the flux calibrators with assumed flux densities of 17.23 Jy at 1280 MHz and 21.03 Jy at 
610 MHz, respectively. The phase calibrators were 2052+365 at 1280 MHz and 1924+334 at 610 
MHz observations, with a bootstrapped flux density of 9.59 Jy and 6.03 Jy, respectively. The data 
analysis was done using AIPS. The estimated uncertainty of the flux calibration is within 7% at 
both frequencies. The image of the field was formed by Fourier inversion and cleaning algorithm 
task IMAGR. The high resolution images (~ 3".3 x 3".0 at 1280 MHz and ~ 5".5 x 4". 8 at 610 
MHz) were made with a Briggs weighting function (robust factor = -1 to 0) halfway between 
uniform and natural weighting. Few iterations of self-calibration were carried out to remove the 
residual effects of atmospheric and ionospheric phase corruptions and to obtain improved maps. 
The system temperature correction was done using sky temperature from the 408 MHz map of 
Haslam et al. (1982). A correction factor equal to the ratio of the system temperature toward the 
source and flux calibrator has been used to scale the deconvolved images. 



3. Other Available Data Sets 

3.1. Near-Infrared Data from 2MAS S 

For the IRSF photometric calibration, study of the bright end of the stellar population and as- 
trometric calibration, we used 2MASS PSC (Cutri et al. 2003) in the J (1.25 yum), H (1.65 jum) and 
(2.17 yum) bands. The 2MASS data have positional accuracy of better than 2". Our source selec- 
tion was based on the 'read-flag' values of "1", "2" or "3", which generally indicate the best quality 
detections, photometry and astrometry ( |http://www.ipac.caltech.edu72mass/reieases/aiisky/doc| ). 

3.2. Mid-infrared data from Spitzer 

The archived MIR data for the region were obtained with the Infrared Array Camera (IRAC; 
Fazio et al. 2004a) on board the Spitzer Space Telescope. The Basic Calibrated Data (BCD) images 
were downloaded from the Spitzer Space Telescope Archive using the Leopard package. The IRAC 
observations of the region were taken on 2005 October 22 (Program ID 20778, AORkey 15154944: 
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A Young Stellar and Protostellar Census of Galactic Ultracompact H ii Regions, PI: S. Carey). Mo- 
saics were built at the native instrument resolution of I". 2 per pixel with the standard BCDs using 
the MOPEX (Mosaicker and Point Source Extractor) software provided by Spitzer Science Center. 
To perform photometry with the DAOPHOT package provided within the IRAF environment, we 
first converted the pixel values of IRAC images, which are in MJy/sr units to DN/s using the conver- 
sion factor given in the "IRAC Data Handbook" version 3.0 (http://ssc.spitzer.caltech.edu/irac/dh/). 
We performed aperture photometry on the IRAC images with a source detection at the 5a level 
above the average local background. Due to the crowded nature of the field, an aperture radius 
of 2 pixels and a sky annulus extending from 2 to 6 pixels were used. The zero point magnitudes 
and appropriate aperture corrections were taken from the the "IRAC Data Handbook" version 3.0 
for the final photometry. The detected sources were examined visually in each band to clean the 
sample of non-stellar objects and false detections around bright stars. We accepted good detections 
as those sources which have photometric uncertainties less than 0.2 mag. 



3.3. Mid-infrared data from MS X 

The Midcourse Space Experimen^iMSX) surveyed the entire Galactic plane within \b\ < 5° 
in four MIR bands centered at 8.28 //m (A), 12.13 jxm (C), 14.65 //m (D) and 21.34 jxm (E) at a 
spatial resolution of 18". 3 (Price et al. 2001) and a global absolute astrometric accuracy of about 
1".9 (Egan et al. 2003). Point sources have been searched around Sh 2-100 region from MSX PSC 
Version 2.3 (Egan et al. 2003). 



3.4. Mid- and far-infrared data from IRAS - HIRES 

The data from the Infrared Astronomical Satellite (IRAS) survey in the two bands (25 and 60 
yum ) for the region around Sh 2-100 were HIRES-proccssed (Aumann et al. 1990) at the Infrared 
Processing and Analysis Center (IPAC), Caltech, to obtain high angular resolution maps. We used 
/7?A5 -HIRES maps to generate the intensity distribution map due to thermal emission from the 
dust component around Sh 2-100. 



"*This research made use of data products from the Midcourse Space Experiment. Processing of the data was 
funded by the Ballistic Missile Defense Organization with additional support from NASA Office of Space Science. 
This research also made use of the NASA/ IPAC Infrared Science Archive, which is operated by the Jet Propulsion 
Laboratory, Caltech, under contract with the NASA. 
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4. General Morphology around Sh 2-100 

Fig. 3 represents a false JHK^ color-composite image (7, blue; H, green; and K^, red) of 
Sh 2-100 region, which reveals prominent nebulosities associated with the marked Hii regions 
except K3-50C1, indicating that the extinction in the molecular cloud towards CI is high enough 
to attenuate the emission at 2.14 jum. A dark lane is also observed in the direction of K3-50C2, 
bisecting its diffuse nebulosity into north and south components, which extend towards K3-50C1. 
A bridge between K3-50D and K3-50A is also seen in Fig. 3, which is not visible in POSS II 
image (Fig. 1), suggesting that the extinction is high towards north-western side of K3-50D. Fig. 
3 also displays several red sources indicating the presence of young stellar sources still deeply 
embedded in the molecular cloud, whereas bluest stars are likely to be evolved/foreground objects 
of the field. The background objects would also appear red, but the locations of high column 
density regions should minimize the contamination of such objects. Fig. 4 shows contours of the 
610 MHz radio free-free emission overlaid onto the IRSF A'^-band image. In the ^g-band image, 
the regions of recent star formation are traced by the bright nebulosity, whereas the H n regions 
detected in centimeter surveys of the cloud trace high-mass star formation. It is well established 
that formation of high mass stars in young complexes is accompanied with the simultaneous or 
subsequent formation of less massive sources, which can be studied in NIR (Persi et al. 2000). 
Our high angular resolution radio map at 1280 MHz is shown in Fig. 5, which was made with 
the robust parameter equal to -1, to see more compact emission features. In Fig. 5 one can notice 
a non-uniform distribution of free-free emission, indicating that Sh 2-100 is morphologically a 
complex region. The low brightness components K3-50E and K3-50F are resolved out in our high 
resolution map at 1280 MHz (Fig. 5). We therefore made a low resolution map (~ 44" x 27") 
using uv range up to 10 kA with a robust parameter equal to 0, to recover the weak diffuse emission 
associated to these sources as seen in 610 MHz map. The resulting map is discussed in §6.5. On 
the other hand, 1280 MHz high resolution map partially resolves two UCH ii regions of K3-50C, 
which are not resolved in our 610 MHz map. 

5. Infrared and Radio Properties of the Sh 2-100 star-forming region 

In order to examine the nature of the stellar populations of Sh 2-100 star-forming region, we 
have used NIR CC and CM diagrams for all the sources detected in JHK^ bands, which are shown 
in Fig. 6. 



5.1. NIR Color-Magnitude Diagram 



In {H - K) versus K CM diagram the nearly vertical solid lines represent zero-age-main- 
sequence (ZAMS) at a distance of 8.7 kpc reddened by Ay = 0, 10, 20, and 30 mag, respectively. 
The parallel slanting lines represent the reddening vectors to the corresponding spectral type. Here, 
we would like to mention that uncertainty in spectral type determination is mainly due to error in 
the distance determination as well as the ^s-band excess of the individual sources. Fig. 6a reveals 
a large nimiber of ZAMS stars earUer than B2 which cannot be explained on the basis of integrated 
Lyman continuum (Lyc) photons from our low resolution radio map. Here, we have incorporated 
low resolution map, because the high resolution (~ 3". 3 x 3".0) observations at 1280 MHz reveal 
only compact sources (see Fig. 5), which will underestimate the total radio flux from the entire 
(~ 8' X 8') region. Since Sh 2-100 is located in the Galactic plane (/ « 70°. 2, I? « r.6) at a 
distance of ~ 8.7 kpc, it is natural that the region could be contaminated by Galactic field main- 
sequence (MS) and giant populations. Therefore in Fig. 6a, we overplotted the locus of giants 
(thick lines at Ay = and 3.3 mag) and supergiants (dashed Une at Ay = 0) using intrinsic colors 
and absolute magnitudes from Cox (2000) for GO to M5 and 09 to M5 stars, respectively. From 
Fig. 6a, one can notice a group of sources located within the reddening strips between Ay = 
and 3.3 mag for giants. Since the region is quite young (see §1), most of these sources must be 
field stars. However, as the region shows strong differential reddening (Roelfsema et al. 1988), it 
is worthwhile to mention that reddened early type MS stars may also fall within these reddening 
strips. The sample may also be contaminated by oxygen-rich asymptotic giant branch (AGB) stars 
and carbon-rich giants, as both these classes of objects can display NIR colors indistinguishable 
from those of normal, reddened stars (Bessell & Brett 1988). To estimate the probable massive 
stellar population, we selected a relatively controlled region (marked as a rectangular box in Fig. 
3) which does not show nebulosity in i?-band (see Fig. 1) as well as in ^s-band (see Fig. 3). The 
region is also devoid of radio free-free emission above 5 mJy at 1415 MHz observation by Israel 
(1976). The 5 mJy flux at a distance of 8.7 kpc for an electron temperature of 10"^ K corresponds 
to the number of Lyc photons ~ 2.9 x 10"*^ s \ which corresponds to a ZAMS spectral type close 
to BO (Panagia 1973). Therefore, if any young massive stars are present in the control region, 
their ZAMS spectral types must be later than BO. The sources detected in the control region are 
shown in Fig. 6b and their positions contradict the above fact, hence we believe that the stellar 
population in our studied region is significantly contaminated by field stars. In comparison to the 
control region, a significant population in the Sh 2-100 region shows red H - K colors, possibly 
the young population of the region. 
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5.2. NIR Color-Color Diagram 

To further characterize the red sources, we use J - H/H - K CC diagram as shown in Fig. 
6c. The solid thin and thick dashed curves taken from Bessell & Brett (1988) represent the locus 
of early MS and late giant stars, respectively. The parallel dashed lines are the reddening vectors 
for early MS and late giant stars (drawn from the base and tip of the two branches). The dotted- 
line represents the locus of Classical T Tauri stars (CTTS, Meyer et al. 1997). The reddening 
vectors are plotted using Ay/Ay = 0.265, Ah /Ay = 0.155, Ak/Av = 0.090 for CIT system (Cohen 
et al. 1981). Since infrared excesses are interpreted as reflecting an early evolutionary stage 
for young stellar objects (YSOs), we classified the CC diagram into three zones (F, T and P) to 
study the nature of sources (for details see Ojha et al. 2004 b,c). The sources in the "F" region 
are generally considered to be either field stars. Class III objects or Class II objects with small 
NIR excess (so-called weak-line T Tauri stars). "T" sources are considered to be mostly Class II 
objects with large NIR excess (CTTS) or reddened early type ZAMS stars, having excess emission 
in ^s-band . There may be an overlap in NIR colors of the upper-end band of Herbig Ae/Be 
stars and the lower-end band of T Tauri stars in the "T" region (Hillenbrand et al. 1992). The 
sources in "P" region are most likely Class I objects (protostar-like objects), having circumstellar 
envelopes. A significant overlap between protostars and T Tauri stars is also seen in CC plane 
by Robitaille et al. (2006). The Galactic evolved stars, star-forming galaxies and active galactic 
nuclei (AGNs) display similar colors as those of YSOs. The thermal emissions from heated dust 
can also be identified as a point source and may appear to the right of the "F" region like YSOs. 
Therefore, NIR CC diagram alone is not an efficient tool to distinguish YSOs from other variety 
of objects, which requires spectroscopic information of individual objects. Due to non-availability 
of spectroscopic observations, we used a statistical approach to identify the probable YSOs in the 
region by comparing the distribution of sources in NIR CC diagram of the Sh 2-100 region with that 
of the control region as shown in Fig 6d. Comparison reveals that sources lying in the "T" region 
and redward of "T" region" with (J-H) > 0.8 in Fig 6c are likely to be the young stellar population 
of the Sh 2-100 region. We do not discard the possibility of a few sources with colors below CTTS 
locus having been missed out in this selection criteria. In this case, we may be underestimating the 
number of YSO candidates. 



5.3. IRAC Color-Color Diagram 

The circumstellar disk emission dominates at longer wavelengths, where the spectral energy 
distribution (SED) significantly deviates from the pure photospheric emission. Hence, the IRAC 
bands can be used to identify IR-excess emission sources (cf. Allen et al. 2004, Megeath et al. 
2004). However, contamination from extragalactic sources such as poly cyclic aromatic hydrocar- 
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bons (PAH)-rich star-forming galaxies and AGNs, limit the use of this method as both types of 
sources display colors similar to that of YSOs. Various authors have investigated the contamina- 
tion of AGNs and galaxies in the IRAC CC diagrams. We followed the approach of Gutermuth 
et al. (2008) and Stern et al. (2005) to identify galaxies and AGNs in IRAC CC diagrams. The 
zones identified by Gutermuth et al. (2008) in [3.6]-[5.8]/[4.5]-[5.8] and [4.5]-[5.8]/[5.8]-[8.0] CC 
diagrams for PAH-rich star-forming galaxies are marked by solid lines in Figs, la and lb, where 
small dots represent the sources identified in all the four IRAC bands within a ~ 8'x 8' region, 
covering the FOV of our NIR observations. The sources marked as crosses in Figs, la and lb 
are star-forming galaxies. Stem et al. (2005) found broad-Une AGNs in a specific zone in the 
[3.6]-[4.5]/[5.8]-[8.0] CC diagram using optical spectroscopic observations, however this zone is 
also occupied by protostars (Allen et al. 2004). Similarly, in a nearby star-forming region, Fazio et 
al. (2004b) found that 50% of the sources with [3.6] > 14.5 mag are extragalactic in nature. Since 
our star-forming region is at a farther distance, some of the YSOs in our sample may appear fainter 
than 14.5 mag in the IRAC 3.6 //m band. Our search of AGNs by combining both the methods as 
mentioned above, results no contamination of AGNs in our sample of YSOs. After identifying the 
extragalactic contamination we used IRAC [5.8]-[8.0]/[3.6]-[4.5] CC diagram, shown in Fig. 7c 
to identify the IR-excess sources (i.e., YSOs). In the CC diagram, the regions occupied by Class 
II, Class I/II and Class O/I objects are indicated following the classification scheme by Allen et al. 
(2004). However, as noted by Whitney et al. (2003b), the SEDs of highly inclined Class II sources 
closely resemble the SEDs, and therefore, the colors of Class I sources (Whitney et al. 2003b; 
Allen et al. 2004). These facts make the boundary between Class I and Class II sources some- 
what indistinct and also Class II objects hidden behind large quantity of dust can also be shifted 
to Class I zone. Here, all the sources with [5.8]-[8.0] > 0.4 and [3.6]-[4.5] > 0.1 are assumed to 
be probable YSOs after eliminating star-forming galaxies. On the basis of the above criteria, the 
[5.8]-[8.0]/[3.6]-[4.5] CC diagram yields 27 YSOs, which are marked as squares in Fig. 7c. The 
search of YSOs on the basis of [5.8]-[8.0]/[3.6]-[4.5] CC diagram could be aff'ected by the lower 
sensitivity of IRAC 5.8 and 8.0 jim bands to detect faint sources. As 4.5 ixm emission traces hot 
circumstellar dust and is less aff'ected by RA.Hs, we used a combination of H, K, and IRAC 4.5 
fim bands to detect IR-excess sources. Fig. Id shows the H - K/K - [4.5] CC diagram. The solid 
curve in the diagram represents the locus of MS late dwarfs (Patten et al. 2006) and the dashed line 
denotes the reddening vector from the tip of M6 dwarf, using the average extinction law by Fla- 
herty et al. (2007). The stars with IR-excess are located to the right of the reddening vector and are 
marked as circles. The squares represent the YSOs selected from the IRAC [5.8]-[8.0]/[3.6]-[4.5] 
CC diagram. Since most of the field stars in the K/H - K CM diagram (see Fig. 6b) have colors 
H-K < 0.35, we therefore adopted a conservative approach to select IR-excess sources in Fig. Id 
by considering sources having H - K > 0.35 and which are also lie to the right of the reddening 
vector. As can be seen in Fig. Id the majority of the YSOs selected on the basis of Fig. 7c fall 
above these color cutoffs. 
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Combining IRAC photometry with NIR observations, we find a total of 150 likely YSO can- 
didates with circumstellar disks, that include Class 1 and Class 11 sources, which indicate a site of 
active star formation. The spatial distributions of these YSO candidates are used to study the star 
formation scenario in this complex, along with their correlation to the associated gas and dust (see 
§7.2). 



5.4. Ionized gas distribution 

There are several radio interferometric studies, mainly towards the regions A, B, C and D 
in the frequency range of 5-15 GHz (Harris 1975; Kurtz et al. 1994; DePree et al. 1994) with 
resolutions better than 3". 5, available in the literature. Our radio continuum maps at 610 MHz 
& 1280 MHz are shown in Figs. 4 and 5, with a resolution of 5". 5 x 4". 8 and 3".3 x 3".0, re- 
spectively. To our knowledge these observations are the highest angular resolution images at low 
frequency to date. The 1280 MHz map shows the radio emission from the five components A, B, 
CI, C2 and D. The components CI and C2 are embedded in an envelope of more diffuse emission, 
whereas the component A shows extended emission in northwest and southwest directions. The 
extended emission in the southwest appears to connect to the low brightness emission seen to the 
north of K3-50D. The measured peak position, integrated flux density and angular diameter, de- 
convolved from the beam for roughly spherical compact sources A, C2 and CI, were determined 
by fitting a two-dimensional elliptical Gaussian using the task JMFIT in AlPS. For irregular, dif- 
fuse extended sources D and B, these observed parameters were obtained using the task TVSTAT 
in AlPS. The flux densities and dimensions for these extended sources were obtained by integrat- 
ing over the approximate angular size of the contour map above 3cr (where cr is the rms noise in 
the map). These observed estimations are given in Table 3. An error of 10-15 % is expected to 
be associated with these values due to the combination of the rms noise in the map (cr), the flux 
calibration error (within 7%) associated with source flux (S) and the uncertainty in flux extraction 
of the source (cTgource)- The combined error is estimated as; cts = sjo-^ + (0.07 Sy + cr^j^j.^^. Here, 
we expected a negligible amount of the missing flux from the compact sources, however in the case 
of extended sources, the uncertainty in the flux estimation could be high, hence the quoted error 
for such sources should be considered as a conservative lower limit. The integrated flux and size 
of the cores of C 1 and C2 should be considered as an upper limit, since the cores are embedded 
in an envelope of difliise emission and the present angular resolution is insufficient to resolve CI 
and C2. It is not known whether the extended emissions of the envelopes are either ionized by 
the sources in the cores and/or (partly) by the nearby external source. The radial velocity of the 
extended emissions and the compact cores are required to know their physical association. The 
images show that all the sources have complicated morphology. The determination of important 
physical parameters is hampered due to the complicated morphology of the sources because these 
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parameters, such as density and surface brightness along the line of sight, depend on the source 
structure. In the present study, we calculate the physical parameters assuming the simplest geome- 
try of uniform spherical model. The physical parameters are such as: He, the rms electron density; 
EM, the emission measure; Tc, the optical depth; Niyc, the number of Lyman continuum photons 
for a spherically symmetric H u region; and the spectral type of the ionizing source, assuming a 
single ZAMS star from the N^yc. These parameters were derived from the measured angular di- 
ameter and integrated flux density using the formulae given by Panagia & Wahnsley (1978) and 
Mezger & Henderson (1967) for spherical, homogeneous nebulae. We have assumed an electron 
temperature of 10"^ K and source size as the geometrical average of the major and minor axes, as 
Usted in Table 3, to derive the physical parameters. The derived physical parameters are Usted 
in Table 4. However, the uncertainties are associated with the values given in Table 4 because 
of our limited knowledge of the source structures, dust absorption of Lyc photons (Garay et al. 
1993) at low frequency for CH n and UCH ii regions (e.g., K3-50A & K3-50C) and missing flux 
from the low density components (e.g., K3-50D & K3-50B). As the gas seems to be clumpy in 
nature, the rms electron density of the H n regions is expected to be lower than the density at the 
peak positions. The sources K3-50C1 and K3-50A are optically thick at 1280 MHz, therefore the 
spectral type, electron density derived by this method should be considered as the lower limit. The 
spectral types of the ionizing sources for the components K3-50B and K3-50D may be earlier than 
the present estimation possibly due to missing flux from the low density components. The spectral 
types suggest that all the sources are ionized by massive late 0-type ZAMS stars and their physical 
parameters point to a different evolutionary class of H n regions (Kurtz et al. 2000). 



5.5. Candidates for the ionizing sources of the H n regions 

The OB type stars associated with K3-50 region are significantly contaminated by the lumi- 
nous field stars (see §5.1). Therefore, in order to identify the probable ionizing candidates of the 
individual H ii regions, we used our JHK^ catalogue. We searched for ionizing sources in a cir- 
cular area of radius ~ 0.5 pc from their radio-emission peak, except for the source K3-50D which 
shows core-halo morphology (see §6.1). We identified the ionizing source of K3-50D using opti- 
cal spectroscopy and discussed in section 6.1. The probable sources are shown in (7 - H) versus 
{H - K) CC diagram (Fig. 8a) as well as in {J - H) versus J CM diagram (Fig. %b). Fig. 9 
shows ^s-band image with the sources located within the radius ~ 0.5 pc of each H n region, along 
with the ionizing source of K3-50D. The sources are labeled with numbers corresponding to the 
H n regions. In Fig. 8fl, sources (Al, A2, C2 and F2) show NIR excess. Therefore, the spectral 
type estimation on the basis of (7 - H) versus 7 CM diagram should be considered as an upper 
Umit for these sources. In order to distinguish the Galactic field stars from the OB stars, we further 
constrain our analysis by using a reddening-free quantity Q = (J - H) - 1.70 (// - )> which is 
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~ for early type stars with no NIR excess and ~ 0.4-0.5 for Galactic field stars and negative for 
NIR excess sources (Comeron & Pasquali 2005). As the sources (Al, A2, C2 and F2) show con- 
siderable NIR excess and may be the candidates for ionizing sources, we selected the sources with 
Q value < 0.3 as possible ionizing sources. It is worthwhile to point out that the sample identified 
using the criterian Q < 0.3 may also be contaminated by objects such as; AGB stars, carbon stars, 
foreground F-type stars and background A-type giants (Comeron & Pasquali 2005). Since we are 
looking for early type stars within the radio emitting region in a circular area of radius ~ 0.5 pc 
around the radio emission peak, we presume the probability of contamination is less, however this 
can not be ruled out. Moreover, in(J-H) versus J CM diagram, we selected luminous sources that 
have spectral type earlier than B3, which further reduces the possibility of contamination. The es- 
timated photometric spectral types for the probable ionizing sources (details are discussed in their 
individual sections) with Q < 0.3 along with the spectral type estimated from radio observations 
at frequencies in which these H ii regions are optically thin, are mentioned in Table 5. Here, it is 
noted that the ZAMS spectral types given in Table 5 are taken from Panagia (1973), corresponding 
to the Niyc photons reported by DePree et al. (1994) from 14.7 GHz observations. It is evident 
from Table 5 that all the H ii regions are excited by early-B to late-O type of stars, which confirms 
the estimation by radio observations, except the source A (see §6.2 for further discussion on source 
A). 



6. Nature of Individual Sources 
6.1. K3-50D 

The morphological features of the nebula at different optical and NIR bands are shown in 
Fig. 10. The ionized gas is seen in Hor with an extended ring, whereas [S ii] image shows a limb 
brightened shell towards the western side. The Bry image traces a semi-circular shell of ionized 
gas with the shell opening towards the north-western side. The same morphology is also seen in 
IRSF i^s-band image. In the IRSF A^s-band image, two bright sources are visible within the ring of 
ionized gas. The star situated at the center of nebula is the main ionizing source of ZAMS spectral 
type ~ 05 (§5.5), whereas the colors (J - H ^ 0.32, H - K ^ 0.20) of the second source situated 
at 12" towards the south-western direction from the central star suggest a reddened giant, which is 
supported by its Q value of ~ 0.38. A plausible model for this type of morphology could be the 
displacement and ionization of the gas due to the stellar wind produced by the massive star at the 

o 

center of the nebula. The one-dimensional spectrum of the star in the range 3950-4750 A is shown 
in Fig. lla. We classified the spectrum by comparing the observed spectrum with those given 
by Walbom & Fitzpatrick (1990). Besides hydrogen Balmer lines (H^, H,5, H^,, etc), the spectrum 
displays prominent He ii line at X 4200 A, and strong He n line at X 4686 A. These features indicate 
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that the star should be an O-type of luminosity class V. In the case of O-type stars the classification 

o o 

of different subclasses are based on the strengths of the He ii lines at A 4541 A and A 4471 A. 
The equal strength of the lines implies a spectral type of 07, whereas the greater the strength of 

o o 

A 4541 A He ii line as compared to A 4471 A line implies a spectral type earlier than 07. If the 
strength of He i + He u A 4026 A < A 4200 A He u , it favors classification earlier than ~ 05.5. 
Comparing the strengths of these lines we adopt the classification as ~ 04V, which is in agreement 
with the estimation by Georgelin (1975). Roelfsema et al. (1988) estimated a visual extinction of 
2 mag towards K3-50D from the ratio of expected Bror emission predicted from radio emissions 
to the observed Bror emission, whereas Persson & Frogel (1974) estimated a visual extinction of 
~ 4.4 mag using the spectrophotometric measurements of hydrogen lines of the nebula associated 
with the K3-50D. To estimate the extinction towards the region, we use the (U - B)/(B - V) CC 
diagram. Morphologically K3-50D shows a bright star at the center and a partial ring of ionized 
gas of diameter ~ 40", therefore we restrict our optical analysis within a circular area of diameter 
~ 50" to avoid field star contamination. The (U - B) versus (B - V) CC diagram of the region 
is shown in Fig. lib. To match the observations, the ZAMS by Schmidt-Kaler (1982) is shifted 
along the reddening vector having a normal slope of |^|r|| = 0.72. The CC diagram indicates a 
minimum reddening of E(B- V)mm = 1-36 mag for the stars within ~ 50" region. Using the relation 
Ay = 3.1 X E(B - V), the estimated value of Ay comes out to be » 4.2 mag. We also estimated 
the average visual extinction for the region, using stars with spectral type earlier than AO by 'Q' 
method (Johnson & Morgan 1953). The intrinsic colors of these stars were obtained from the 
relation (B - V)o = 0.332 x Q (Hillenbrand et al. 1993). To check the bias of the contamination of 
other field sources to this value, we estimated the visual extinction from the ionizing source of K3- 
50D. The O-type MS star shows a narrow range in their intrinsic colors. Adopting an intrinsic color 
(B - V)o = -0.33 (Cox 2000) for the ionizing star of K3-50D and using the measured photometric 
color (B-V) = 1 .04, we derived an excess E(B-V) = 1 .37, which in turn gives Ay w 4.2 mag, which 
agrees well with the estimation based on the CC diagram. Since the distance estimates to the region 
are uncertain and vary from 4.5 to 12.5 kpc (Israel 1976), the present spectral classification of the 
ionizing star and reliable estimation of extinction together with the present apparent magnitude of 
the ionizing source enables us to derive its distance. Using the present spectral classification (04V) 
along with Ay « 4.2 mag, V = 13.31 and My -spectral type calibration table of Vacca et al. (1996), 
we estimated a distance of ~ 8.5 kpc for the K3-50D region. However, the spectroscopic distance 
determination is strongly dependent on the assumed value of My and the values of My in case of 
04V star vary from -5.4 (Schmidt-Kaler 1982) to -5.9 (Conti 1988) in the literattire, which in turn 
can vary distance from 7.9 kpc to 10 kpc. However, recent investigation by Martins et al. (2005) 
found considerable shift in calibration of effective temperature in early type stars. The My value 
of 04V star (My = -5.5) by Martins et al. (2005) results in a distance of 8.6 kpc to the region. 
These estimations are consistent with the average kinematic distance (~ 8.7 kpc) towards K3-50 
and supports its association with the W58 molecular cloud. 



-17- 



6.2. K3-50A 

K3-50A is a CH ii region, ionized by a single 05.5V star and its angular size is ~ 4".l x 3".4 
at 14.7 GHz (Kurtz et al. 1994). The radio continuum peak is coincident with 10 jum peak, 
whereas the associated optical emission (see §1) is located to the south-west (~ 2.5") of the radio 
peak (Wynn- Williams et al. 1977) and the region is possibly ionized by multiple sources (Okamoto 
et al. 2003). We detected two NIR sources in all the three bands within ~ 0.5 pc radius of the radio 
peak. The source Al shows a large IR-excess (H-K » 3.35) and is also more luminous (K » 9.1 1) 
than any reddened 0-type stars of the complex. The second bright source A2 is located outside the 
radio size (~ 4") of K3-50A, hence the probability of this source as one of the ionizing candidates 
is rather low. Therefore, Al is the most probable ionizing source of K3-50A. We speculate that 
all the point sources close to the nebula (e.g., A2) may be strongly contaminated by the intense 
/Ts-band nebulosity of K3-50A. The crowding in the region has already been noticed by Okamoto 
et al. (2003), therefore, the spectral classification on the basis of NIR photometry may be affiected 
by source contamination along with bright nebulosity and will yield an earlier spectral type. 



6.2.1. Spectrum of the Nebula 

Fig. 10(2 shows the continuum- subtracted Ha -I- [N ii] image of the optical nebula associated 
with K3-50A. The spectrum of the nebula is presented in Fig. 12. The nebula shows forbidden 
emission lines as seen in the case of H n regions. We obtained the emission line fluxes by Gaus- 
sian fitting of the line profiles and by direct integration of the fluxes over the local continuum 
with SPLOT routine of the IRAF package. Uncertainties in these line intensities depend upon the 
strength of the lines and blending with adjacent Unes. The measurement errors for the isolated 
bright lines (fluxes larger than ~ 10"^"^ ergs cm"-^ s~^) are expected to be less than 10%, whereas 
for the weaker and blended Unes like Ha + [N ii] and [S ii], the measurement uncertainties can 
be as large as 20%. Our measurements are consistent with that of Persson & Frogel (1974) for the 
common lines. 

The physical parameters of the nebula can be derived from the diagnostic of line ratios with 
certain assumptions. Since our spectrum contains emission lines of hydrogen, their relative inten- 
sities can be used to estimate the reddening for the region. We obtained the reddening coefficient 
CiHjs) « 3.16, by comparing the observed ratio of Ho'/H/j with the theoretical ones calculated by 
Osterbrock (1989), for the temperature lO'* K and density 10"* cm~^, using Case-B recombination 
theory and extinction law of Cardelli, Clayton & Mathis (1989, hereafter CCM 89). Here, we 
would like to mention that the presence of extinction increases with the observed ratio, because 
extinction affects the H^ wavelength more than the Ha wavelength and also the ratio depends on the 
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geometry of the dust and gas cloud. However, Persson & Frogel (1974) found variable extinction 
within the region by comparing the radio, infrared and optical emission fluxes. Wynn-WilUams et 
al. (1977) found that these emissions are not coincident and the location of peak intensity moves 
from the radio center towards the southwest direction as the wavelength decreases. The position 
of the optically visible nebula is approximately 2". 5 away from the peak of radio position. On the 
basis of a deep 9.7 yum absorption feature in the MIR spectra, it is interpreted that K3-50A is partly 
obscured by a dense dust cloud (Thronson & Harper 1979 and references therein), which sup- 
ports similar conclusion of Wynn- Williams et al. (1977). To explain the strong extinction gradient 
across the region, Wynn- Williams et al. (1977) proposed a possible geometry that the K3-50A is a 
partially ionization bounded CH n region formed near the edge of an obscuring cloud and is now 
eating its way to the less dense region at the edge, like the Orion Nebula H ii region. 

Assuming that the optical lines represent the physical properties of local environment, we re- 
tain our analysis with Ay ~ 6.79 mag, which has been derived using the extinction law Ay ~ 2. 145 
X C(H(^) (CCM 89). All the line fluxes were dereddened using Ay ~ 6.79 mag and extinction law 
of CCM 89. Table 6 lists the observed and reddening corrected emission line intensities relative to 
H^. We estimate the electron density ~ 0.67 x 10"^ cm"^ using 5-level atom program by Shaw & 
Dufour (1994) from the [S n] (A 6111 / A 6731) line ratios, assuming Te = lO'* K. This assumption 
is not critical because the density calculated from [S u] lines are nearly independent of Tg. The 
[S ii] line ratio (~ 0.54) lies close to the high density limit, probably due to the coUisional quench- 
ing of forbidden Unes, hence the estimated electron density should be considered as a conservative 
lower limit though it might be higher. It is noted that the density derived from the [S n] lines is 
representative of a nebula ~ 2".5 away from the high-density core, traced by high frequency high 
resolution radio observations. Therefore, the electron density estimated with high frequency radio 
observations is expected to be higher than density from [S ii] lines. We estimate the electron tem- 
perature ([N ii] ) from the temperature sensitive intrinsic [N ii] (6548 A+ 6583 A/ 5755 A) line 
intensity ratio, which turns out to be ~ 10496 K. However, it should be noted that 5755 A line is 
quite faint as seen in the spectrum. A roughly linear, increasing trend in [O III] (4959 A+ 5007 
AyHj} with stellar temperature is expected, because the nebular volume containing doubly ionized 
oxygen increases with the higher ionizing flux emitted by early spectral type stars. By comparing 
the [O III] (4959 A+ 5007 A)///^ value ~ 1.79 of K3-50A, with that of values in the range 1.71 to 
1.86 for the Galactic H n regions ionized by a single star (Kennicutt et al. 2000), a rough estima- 
tion of stellar temperature of the ionizing star hotter than 35900 K and cooler than 42700 K can be 
obtained. Thus, from the stellar temperature, a MS spectral type between 09 to 06.5 (Vacca et al. 
1996) is expected as the probable exciting source of the nebula. 
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6.2.2. Sub-millimeter emission from cold dust 

As star-forming regions are optically thin for wavelengths longer than ~ 200 //m, thermal 
dust emission can be used to determine the physical parameters such as mass and density of star- 
forming regions. Continuum sub-millimeter fluxes at 850 jxm and 450 jxm for the K3-50A were 
obtained from Thompson et al. (2006). The integrated flux densities at 450 iira and 850 jira are 
256 + 79.3 Jy and 37 + 4.3 Jy, respectively. The peak flux densities are 53.6 + 16.1 Jy/beam and 
13.6 + 1.4 Jy/beam at 450 and 850 fxra, respectively. The 850 jira map is shown in Fig. 13. The 
bright continuum emission is concentrated towards the central region of K3-50A, showing a core 
peaking at 0^2000 = 20''01'"45^, ^2000 = +33°32'44", with extended emission along the northeast 
and southwest directions. For optically thin emission, the dust mass can be estimated from the 
standard method of Hildebrand (1983) and the total mass (dust plus gas) of the cloud is given by 

M=^, (1) 

Where Fy is the flux density at frequency v, By{T^ is the Planck function evaluated at frequency 
V and dust temperature T^. The parameter Cy is a mass conversion factor combining both the 
dust-to-gas ratio and the frequency-dependent dust opacity Ky and d is the distance to the source 
in kpc. Various values for Cy are quoted in the literature; here we have adopted a value for Cy 
of 50 g cm~^ at 850 //m, following the method of Kerton et al. (2001), which is appropriate for 
dense molecular cores. The integrated flux over the region at 850 jum gives the mass of the cloud 
as ~ 7800 Mq for an assumed average dust temperature of 30 K for the entire cloud. The main 
error in this estimate is likely due to the uncertainty in the dust temperature. Because of non-linear 
dependence of the Planck function on temperature, a small change in temperature leads to a large 
change in the calculated dust mass (e.g., a 5 K error in temperature leads to an error of ~ 40% 
in mass estimation). Hence we consider the present mass estimate from the dust emission to be 
accurate only within an order of magnitude. 



6.2.3. Spectral Energy Distribution 

The star Al shows large infrared excess (Fig. 8a) and is located close to the peak of the molec- 
ular core having weak bipolar outflow (Phillips & Mampaso 1991), a phenomenon that is thought 
to be closely related to the accretion process. Molecular outflows are common phenomenon in the 
early stage of massive star formation prior to formation of UCH n region (e.g., Beuther et al. 2002) 

and also during the UCH 11 phase (e.g.. Shepherd & Churchwell 1996). To investigate the nature of 
Al, we constructed SED using the recently available grid of models and fitting tools of Robitaille 
et al. (2006, 2007). The models are computed using a Monte-Carlo based radiation transfer code 
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(Whitney et al. 2003a, b), using several combinations of central star, disk, in-falling envelope and 

bipolar cavity, for a reasonably large parameter space. To produce the observed SED, we have used 
the data from the following: 1.25-2.17 //m (Alvarez et al. 2004), 3.5 and 10.1 /im (Neugebaure & 
Garmire 1970); 8-21 /urn, MSX PSC (Egan et al. 2003); 450 and 850 /urn (Thompson et al. 2006). 
The SED is shown in Fig. 14. Here, we would like to mention that the angular resolutions for 
most of the data are not sufficient to derive precise parameters for a single YSO. We set MSX and 
sub-mm fluxes as upper limits to account for possible contribution from other sources that might 
be included within the large beam of the observations. We further set 20 to 10% error in NIR 
and MIR flux estimates due to possible uncertainties in the calibration, the extinction and intrinsic 
object variability. A distance range of 7 to 10 kpc and the visual extinction in the range from 
15 to 30 mag are used as input parameters to the model to construct the SEDs. The parameters 
obtained from the weighted mean and standard deviation of all the models with - Xhest^ < 3 
(per data point), weighted by inverse square of^^ of each model, suggest the associated source is 
a massive YSO with mass = 42 + 4.4 Mq and of total luminosity Log(L^) = 5.6 + 0.2 Lq. The 
models also predict that the source is deeply embedded behind 17.5 + 1.8 mag of visual extinction 
and has an envelope accretion rate of Log(Me„„) = -3.8+ 0.2 Mq yr~^ with no disk. The SEDs of 
high mass protostars are dominated by the envelope flux, where the disk flux is almost embed- 
ded within the envelope flux. Therefore, the models with envelopes alone and the models with 
disks and envelopes successfully fit the same SED (Grave & Kumar 2009). It is to be noted that 
these models are valid for isolated objects. Therefore, considering the presence of multiple point 
sources claimed by Okamoto et al. (2003), the above parameters should be taken with caution. 
Nevertheless these models provide quantitative information about the nature of the object. 



6.3. K3-50B 

The ^s-band image of K3-50B shows unusual distribution of dust within the H n region with 
an elongated hole from east to west (see Fig. 3). The morphology of 610 MHz map (Fig. 4) 
suggests a blister type H n region, where a sharp ionization bounded edge can be seen towards the 
western direction. Towards the east, the gradient of flux density indicates a density-bounded side 
of the blister morphology. The high resolution image at 1280 MHz (see Fig. 5) also shows non- 
uniform clumpy distribution of ionized gas. The complex structure of radio emission suggests that 
it may have been produced through excitation by multiple stars (e.g., Ojha et al. 2004a; Garay et al. 
1993) or presence of density inhomogeneities within the H n region excited by a single luminous 
star. We detect four sources (marked in Fig. 9) within a radius of ~ 0.5 pc from the center of the 
nebula. The source B4 is a bright infrared source (K « 11.27, H - K 0.88) with Q « 0.06, 
whereas the source B3 {K « 12.47, H - K = 0.31) has Q » 0.31. The source B3 is bright in 
the optical wavelengths and being situated at the central region of the nebula gives an impression 
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of it being one of the ionizing candidates. The optical spectrum of the source B3 in the range 

o 

3950-4750 A is shown in Fig. 15. Due to low S/N ratio, the exact spectral classification of the 
source is rather difficult, but the absence of absorption lines (He n at A 4200 & A 4686 A, He i at 
A 4026 & A 4471 A, Si n at A 4128 & A 4130 A and Mg u at A 4481 A), suggest that source B3 
is not an OB-type star. The source Bl has Q « 0.49, along with H - K x 0.61, which suggests 
a reddened giant rather than a ZAMS star. The positions of the sources Bl and B3 on the CC 
diagram (Fig. 8a) fall close to the giant locus, confirming the above argument. The source B2 
with 2 ~ 0.12, indicates that it could be one of the ionizing candidates of ZAMS spectral type of 
B2-B3, but its low visual extinction (Ay » 3.8 mag) in comparison to Ay » 26 mag (Roelfsema et 
al. 1988) towards K3-50B suggests a field star. The source B2 seems to have spectral type later 
than B2. As the star of spectral type later than B2V emits < 2 % of its total output energy in terms 
Lyc photons per second, therefore the source B2 cannot be one of the main ionizing candidates of 
K3-50B. The position of the source B4 on NIR CM diagram (Fig. Sb) shows the characteristics 
of ZAMS spectral type ~ 05, which is consistent with spectral type (~ 06-05.5) derived from 
Lyc fluxes from radio observations, therefore B4 can be the most probable ionizing candidate for 
K3-50B. The spectral type estimation of B4 is not biased by excess emission as it does not possess 
NIR excess, therefore the complex structures of the K3-50B region are probably due to density 
inhomogeneities rather than due to multiple early-type sources in the region. 



6.4. K3-50C 

K3-50C is located ~ 2'. 5 northeast from K3-50A and consists of two radio continuum sources 
CI and C2. These sources are separated by ~ 15" and CI is stronger than C2 at 8.4 GHz (Kurtz et 
al. 1994) and 15 GHz (DePree et al. 1994) radio wavelengths. However, below 20 jum C2 is much 
stronger than CI and CI is barely detected at 20 //m and is not seen at shorter wavelengths (Wynn- 
Williams et al. 1977), indicating that CI is more deeply embedded in molecular cloud of high 
visual extinction as compared to C2. C2 is located at the center of a prominent ^g-band nebulosity 
in K3-50C region, where we detect a bright NIR point source close to the radio peak. The positions 
of the source C2 on CC and CM diagrams (Figs. 8a & Sb) indicate a reddened ZAMS spectral type 
between BO - 09.5 with NIR excess. The NIR spectral type estimation agrees well with the ~ 09 
ZAMS spectral type derived from radio continuum observations. The position of CI falls on the 
dark lane that bisects ^s-band diffuse emission into north and south components (Fig. 3). A large 
variation of extinction is found in the direction of CI. At the peak of radio continuum emission 
no Bror was detected, indicating a visual extinction of > 190 mag (Roelfsema et al. 1988), which 
is the possible cause for the non-detection of ionizing source(s) in our infrared survey around the 
radio peak. The Niyc photons s~^ from the region suggest that CI is ionized by a single 08-09 
ZAMS star. The spectral type estimation on the basis of radio continuum photons is an indirect 
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approach in the absence of estimation based on direct stellar radiation from the star(s). Therefore, 
here the possibility that more than one star might be responsible for the ionization of CI cannot be 
excluded, since massive star formation usually proceeds in cluster mode. In this case, the spectral 
type of the most luminous star must be later type than that derived from the single-star hypothesis. 



6.5. K3-50E and K3-50F 

K3-50E and K3-50F are two faint optical nebulosities separated by 2'.2 (Fig. 1). Our low 
resolution (~ 44" x 27") map at 1280 MHz (see Fig. 16) shows continuum emission from ionized 
gas around both the regions, with K3-50F being relatively brighter than K3-50E. In Fig. 4, bright 
^s-band sources can be seen within the radio contours, which are also visible in the optical (K3- 
50E: V = 17.16, B-V = 1.543; K3-50F: V = 14.95, B-V = 1.282). These sources are located 
within the radius of ~ 0.5 pc from the radio peak, and so are most likely the ionizing sources of 
the nebulae. Their visual extinctions estimated using the 'Q' method (§6.1) suggest that they are 
extincted by visual extinction of 5.6 and 4.8 mag, respectively. 

There are no definite estimations of distance for K3-50E and K3-50F in the literature but 
they are assumed as ~ 8.7 kpc (Israel 1976). Hence, to estimate approximate distance we use 
the following relation given by Comeron & Torra (2001) for an ionization bounded optically thin 
H II region with an electron temperature of ~ 10"^ K : 

D{kpc) « 5.35 X — I ^_^Jif.233 (2) 

where Sy is the flux density in Jy at frequency v (in GHz). Hoi is the dereddened //-band magnitude 
of 'i' th ionizing star. Due to the small sizes of both the regions and our selection method to identify 
possible ionizing sources described in section 5.5, we believe that the two optically bright sources 
are most likely the ionizing sources of the H n regions. Using the above approximate relation, we 
obtain distances of ~ 7.6 kpc and ~ 8.4 kpc for K3-50E and K3-50F, respectively. Here, we assume 
that both the regions are ionization bounded and ionized by a single star. In Eq. 2 we have used 
radio fluxes 0.04 Jy for K3-50E and ^ 0.13 Jy for K3-50F) estimated from the low resolution 
map (Fig. 16), the //-band magnitudes (« 12.75 mag for K3-50E and « 12.25 mag for K3-50F) of 
the ionizing stars and extinction at //-band (« 0.87 mag for K3-50E and « 0.74 mag for K3-50F). 
Although the above approximate distance estimates agree with the kinematic distance (7.9 to 9.3 
kpc) reported towards the direction of K3-50 SFR, the validity of the method is constrained by 
the number of assumptions concerning membership, dereddening, dust absorption of Lyc photons, 
leakage of ionizing photons from the H ii regions and uncertainty in the radio flux estimations. 
The more precise distance can be estimated using observed photometric magnitudes and accurate 
spectral type of the ionizing sources as described in section 6.1 for K3-50D. We estimated Lyc 
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photons per second from the integrated fluxes at an assumed distance of 8.7 kpc and that come 
out to be ~ 1.5 X 10"*^ s^^ and ~ 4.8 x lO'^'^ for the K3-50E and K3-50F, respectively. These 
fluxes correspond to ZAMS spectral types of B0.5-B0 and BO-09.5, respectively, consistent with 
the ZAMS spectral types estimated from J - H versus J CM diagram (see Fig. 8& and Table 5). 



7. General discussion 

7.1. Extinction towards individual regions 

The massive stars (M > SMq) do not have pre-main-sequence (PMS) phase and the object 
reaches ZAMS while still embedded in the molecular cloud (Palla & Stabler 1990). We used 
probable ionizing sources (see Table 5) of the individual H u regions to derive their extinction 
values. We traced the ionizing stars back to their intrinsic ZAMS position in (/ - ^ versus / CM 
diagram (Fig. 8&) along the reddening vector. However, the uncertainty in estimation of extinction 
cannot be ruled out due to excess emission in J and H bands. This could be the case for the ionizing 
star Al of K3-50A region, which shows maximum NIR excess among all the YSOs. The visual 
extinction toward regions A, B, C2, D, E and F is estimated to be 31.2, 16.1, 21.2, 4.4, 5.9 and 
5.5 mag, respectively. Here it is to be noted that the extinction estimates for Al and C2 must be 
towards the higher side as these stars have significant NIR excess. For example, in the case of Al, 
the extinction value derived from the SED fitting (§6.2.3) comes out to be 17.5 + 1.8 mag. The 
above discussion along with the average extinction found towards CI (see § 1) indicate that the K3- 
50 SFR shows variable extinction in the range of 4.4 - 97 mag. A comparison of extinction found 
towards K3-50 to the i^CO(l-O) map by Israel (1980) for W58 complex at a resolution of 2'.3 
reveals the following facts: the position of K3-50C coincides with the ^^CO(l-O) peak, whereas 
the regions K3-50D, K3-50E and K3-50F are located towards the edge, and therefore suffer less 
amount of extinction (see Fig. 17). The relatively better resolution (~ 1') 1 mm continuum map 
of Wynn- Williams et al. (1977) shows mainly two peaks (towards K3-50C and K3-50A) of dust 
emission, indicating that K3-50B is less obscured in comparison with both these sources. This 
agrees well with its low extinction value. 

7.2. Spatial distribution of YSOs 

The IR-excess in the case of YSOs can be due to circumstellar disk/envelope around the YSOs 
or weaker contribution from reflected stellar radiation of the dust emission. In any case IR-excess 
represents the association of young sources in a star-forming region. The spatial distribution of 
probable YSOs viz. Class II YSOs (asterisks). Class I YSOs (triangles), possible YSOs detected 
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only in H and bands with H-K>\.i (filled circles), and selected from H - K/K - [4.5] CC 

diagram (open circles) is shown in Fig. 17. Most of the Class I sources and sources with H - K > 
1.8 are mainly concentrated towards the regions A, B and C2. The detection of sources with H - K 
> 1.8 towards A, B and C2 suggests that the extinction is so large that many YSOs may be too 
heavily extincted to be detected in 7-band. Since most of these sources are also associated with 
/Tj-band nebulosity, they possibly have intrinsic NIR excess as well as local extinction. However, 
the IR-excess sources identified from H - K/K - [4.5] CC diagram are spreaded over the entire 
region with less concentration towards the south-eastern direction. The spatial distribution of these 
IR-excess sources is well correlated with the difiiise emission seen in 8 fxm band of Spitzer as well 
as with the intensity of CO emission. An inspection of the MSX PSC yields six sources within 
the investigated region with positive flux values in A (8.3 )Um), C (12.1 jum), D (14.7 yum) and E 
(21.3 yum) of MSX bands. The positions of MSX point sources are shown in Fig. 17 as crosses. 
The spatial distribution of four sources coincides with the ^s-band nebulosity associated to the 
Hu regions. The position of one source (G070. 2990+01. 5762) falls on the north-eastern side 
of the K3-50A and the other source (G070.3164+01.6493) coincides with the intense compact 
^s-band nebulosity seen between K3-50E and K3-50F (see Fig. 3). An IRAS point source (IRAS 
19597+3327 A) is also detected in the proximity of G070. 3 164+0 1.6493 (hereafter IRAS-B, see 
details in appendix). To know the nature of these sources, we use CC plot as suggested by Lumsden 
et al. (2002) who derived criteria for identifying various Galactic plane sources using the MSX 
bands data. Fig. 18 shows the CC plot for the MSX point sources. The demarcations shown in 
Fig. 18 were derived from known sample of objects by Lumsden et al. (2002), where a significant 
overlap between massive YSOs and CH ii regions is seen in their sample having colors FjxIF^ > 
8 and Fu/Fn > 1. Figure 18 shows that most of our sources fall in the overlap region of massive 
YSOs and CH n regions. To study the nature of these sources we constructed SEDs with the help of 
their NIR and IRAC counterparts. For the construction of SEDs, IRAC sources with photometric 
error less than 0.25 have been used. Fig. 19 displays the SEDs for the sources between 1.25 fxm to 
8 fim. Protostars are classified on the basis of their spectral slopes from their SEDs. The general 
consensus is that rising and flat spectrum sources are, on an average, younger and are in a more 
embedded phase of star formation (Allen et al. 2004). The SEDs of massive stars associated to D, 
E and F show almost flatter spectrum or a minor increasing trend, whereas those of IRAS-B and 
C2 show rising trend at longer wavelengths. From SEDs it appears that the source C2 is relatively 
younger than the other sources, neglecting the eff'ect of inclination, stellar temperature, contribution 
from the scattered Ught and extinction, which can cause diversity in the spectral slopes. 



7.3. Age of the massive stars 



The ages of young clusters are typically derived from the post-main-sequence evolutionary 
tracks, for the earliest members, if significant evolution has occurred. The presence of early O- 
type MS star (e.g., 04V in the case of K3-50D) in a star-forming region would indicate that the 
age of the region should be < 3 Myr, as lifetime of an 04 star is of the order of 3 Myr (Meynet 
et al. 1994). Massive OB stars do not have optically visible PMS phase. As Sh 2-100 region 
contains a few optically detected early-B stars and the contraction time scale for such stars (B3V; 
9M© to BOV; ISM©) to reach ZAMS is approximately few x 10^ yr (Bemasconi & Maeder 1996), 
hence the expected ages of massive stars associated to the Sh 2-100 region must be of the order 
of few Myr. To estimate the ages of ionizing sources of K3-50D, K3-50E and K3-50F, we used 
V/B-V CM diagram. Fig. 20 shows intrinsic V/B-V CM diagram, where the stars are dereddened 
individually using the Q method as discussed in section 6.1. We visually fit solar metallicity theo- 
retical isochrone of 1 Myr by Girardi et al. (2002), assuming a distance of 8.7 kpc for the region. 
However, due to the distance uncertainty associated with the individual sources, uncertainty in the 
estimation of extinction based on Q method and the resolution of the isochrones to distinguish the 
ages between 1-3 Myr, the derived age (~ 1 Myr) should be considered as an approximate estima- 
tion. Probable ionizing stars of the remaining regions are embedded and not detected in the optical 
band. Consequently, their age estimation is rather difficult. Most of these stars show NIR excess, 
which is a signature of the presence of an inner disk. As the time scale for destruction of disks 
ranges from less than 10^ yr for an early O star to ~ 5 x 10^ yr for BOV star (Bik et al. 2006 and 
references therein), the ages of the massive stars associated with K3-50A, K3-50B and K3-50C2 
should be less than 1 Myr. 



7.4. Mass Spectrum of Young Stellar Objects 

Considering the presence of CH n & UCH n regions, astronomical masers and intense gas 
and dust, similar to the W3 Main region studied by Ojha et al. (2004b), we assumed PMS age of 
0.5 Myr to 1 Myr for the YSOs to estimate their masses. To minimize the effect of NIR excess 
emission on mass estimation we preferred /-band luminosity rather than H or K band. Fig. 21 
represents {J - H) versus / CM diagram for the YSOs. The solid and dotted curve lines in Fig. 
21 denote the loci of 0.5 Myr and 1 Myr PMS isochrones by Siess et al. (2000). The soUd and 
dotted slanting lines are the reddening vectors for 0.3 and 4 M© stars for the 0.5 Myr and 1 Myr 
isochrones, respectively. The different symbols represent the YSOs selected from NIR and IRAC 
CC diagram and are the same as in Fig. 17. The reddening vectors for the assumed age suggest 
that the majority of the YSOs have masses in the range of O.3M0 - AM^ indicating that stellar 
population in Sh 2-100 is mainly composed of low-mass YSOs. Although, the assumed age for 
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the region is the subject of uncertainty and may not be accurate as compared to the ages estimated 
with the spectroscopic information of the sources in star clusters, the distibution of YSOs with wide 
colors (Fig. 21) probably indicates the combined effect of variable extinction, weak contribution 
of excess emission in J and H bands and/or sources in different evolutionary stages. Therefore, it 
would be very useful to get radial velocity information of the YSOs to confirm their association 
and ages with spectroscopic information. 



7.5. A possible evolutionary status of the H n regions 

The regions K3-50E and K3-50F contain massive BO-BO. 5V MS stars. These early-type 
stars usually form in the dense core of a molecular cloud and can be seen in millimeter, mid-, 
far-infrared and/or as compact radio emissions. The presence of weak extended radio emission 
(Fig. 4) and lack of dust continuum emission at 25 //m and 60 //m (see Fig. 22) indicate that the 
associated H ii regions are probably in an advanced stage of evolution. The presence of relatively 
intense radio emission at 1280 MHz and dust emission at 25 fxra and 60 jum, respectively, indicate 
that K3-50D is younger or of age comparable with K3-50E & K3-50F. All the ionizing sources 
of these H n regions do not show infrared excess and their detection in the optical band suggests 
their evolved nature. The optical V/V-I diagram suggests that the age of the ionizing sources is of 
the order of 1 - 2 Myr. However, non-detection of ionzing source(s) in optical, relatively intense 
-band emission as well as radio emission as compared to K3-50D suggest that K3-50B should 
be younger than K3-50D. The regions K3-50A and K3-50C (consisting of two sources CI and 
C2) are possibly younger than K3-50B, due to their association with intense ^^CO(l-O) emission, 
compact radio sizes and the ionzing sources possess infrared excess (e.g., C2 and A). The detection 
of OH masers (Baudry & Desmurs 2002) and ionized outflows (DePree et al. 1994; Balser et al. 
2001) in both the cores indicate that they are not dynamically evolved objects and suggest that the 
star formation activity around K3-50C and K3-50A should not be more than ~ few x 10^ yr. The 
component CI is hidden behind a visual extinction > 190 mag, therefore invisible in wavelength 
shorter than 20 yum, whereas the ionizing source of K3-50A is detected in NIR and embedded in 
molecular cloud having visual extinction ~ 17.5 mag. K3-50C1 has lower dust temperature than 
K3-50A, possibly indicating that K3-50C1 is more deeply embedded in the molecular cloud as 
compared to K3-50A (Hunter et al. 1997). The above discussion reveals that K3-50C1 is most 
Ukely the youngest source of the region. 
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7.6. Star formation scenario in W58 

W58 is a molecular cloud complex of an angular extent of about 60 arcmin. The ^^CO dis- 
tribution mainly reveals three distinct components (see Fig. 3 of Israel 1980) and the associated 
SFRs to these components are K3-50, Sh 2-99 (S99) and W58G. Fig. 23 shows 8 /im emission 
from MSX A-band along with the distribution of H i emission around the SFR complex (thick 
solid line) taken from Fig. 6 of Israel (1980). The positions of different sources found towards the 
direction of W58 cloud complex are shown as crosses. The most intense northern CO component 
associated with K3-50 SFR peaks close to K3-50C, where most of the star formation activity is cur- 
rently going on. The western CO component is associated with an optically visible H ii region (Sh 
2-99/G70. 15-1- 1.73) situated at a distance of ~ 11' from K3-50A towards the south-western direc- 
tion. The weaker southern CO component is associated with an H ii region (W58G/G69. 92-1- 1.52) 
which is ~ 23' away from K3-50A. In the case of H n region (W58G), a radial velocity of -65 km 
s~^ has been measured, which differs significantly from the radial velocity - 22 to -25 km s~^ found 
in the direction of W58. This indicates that W58G most likely is not related to W58 cloud complex 
(Israel 1980 and references therein), whereas the radial velocity of -22.9 km s~^ for Sh 2-99 (Blair 
et al. 1975) suggests its association with W58. 

With the help of H i emission observed by Read (198 1) with a resolution of 2'. 5 x 4'. 6, Israel 
(1980) pointed out an H i /H n expanding shell of diameter ~ 75 pc and suggested that the current 
star formation activity in Sh 2-99 and K3-50 region might have been induced by the interaction 
of this shell produced by stellar winds from the older generation stars of W58 cloud complex. 

To recognize such a shell, a potential way is to use the R^.H emissions at 8 fim arising from the 
photodissociation region, which occurs at the interface between the molecular and neutral gas and 
thought to be excited by sub-Lyman UV photons from the massive star(s). Consequently, it can be 
used as a tracer of the neutral shell observed in the H 1 21 -cm line emission at the edge of expanding 
H II region (Deharveng et al. 2005). Fig. 23 shows the 8 /urn R\H emission surrounding the bubble 
noticed by Israel (1980) and roughly coincides with the H i emission in terms of discrete patches. 
The absence of 8 fj.m emission in the interior of the H i shell can be interpreted as the destruction 
of R\H molecules by intense UV radiation from the associated massive star(s). The shape and 
large dimension of H i shell suggest that the stars producing stellar wind should be located inside 
the shell. A search of Simbad database results in a WR star (WR 131) of WN7h (hydrogen rich) 
subtype (Figer et al. 1997), whose position is marked in Fig. 23. The hydrogen rich WN7 stars 
are descendants of massive stars with initial masses above 50-60 Mq (e.g., Crowther et al. 1995). 
Indeed, progenitors of such stars are believed to be core-hydrogen burning Of stars with strong 
stellar winds (Crowther et al. 1995). The presence of WR 131 in the H i shell and its approximate 
distance of ~ 9.1 kpc (Georgelin et al. 1988) comparable with the distance estimates of W58 
complex (~ 7.9 to 9.3 kpc), favors a probable association between the WR 131 and W58 complex. 
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The catalog of Reed (1998) yields a few B-type stars which are projected inside the bubble, but their 
distances are not known so as to confirm their association. If indeed these sources are associated 
then the presence of H i shell can be due to the combined energy input from these sources. 

WR 1 3 1 is not located at the center of the H i shell but at an off-center position, which could be 
due to the combination of an inhomogeneous medium, proper motion of the star, anisotropic winds 
and projection eff'ect along the line of sight. If we assume WR 1 3 1 is the source to produce H i shell 
due to its strong stellar winds and since WR phase is short lived (~ 10^ yr) and the presence of 
WR star limits the maximal age ~ 4.5 - 5 Myr of a cluster (Meynet & Maeder 2005), an upper 
limit of age ~ 4 - 4.5 Myr can be assigned to WR 131, considering the estimated dynamical age 
(~ 10^ yr) of the bubble by Israel (1980). If the age of WR 131 is ~ 4 - 4.5 Myr then one would 
expect the ages of the second generation stars, (if any), at the periphery of the bubble, should be 
less than that of a WR star. The presence of star formation activity, which might have started < 
2 Myr ago is evident at the north-eastern edge of the shell, with the existence of Class I sources, 
astronomical masers, UCH n and CH n regions containing young OB stars of ages 1-2 Myr. The 
presence of seven H ii regions at different spatial locations and possibly at different evolutionary 
stages suggests star formation activity has proceeded at multiple sites in the cloud. The presence 
of WR 131 star of age ~ 4 - 4.5 Myr and young massive stars still in accretion phase (e.g., CI 
and A), associated with the same complex, indicates distribution of stars of different ages in the 
region. Israel (1980) suggested that at present the shell is colliding with a larger H i cloud to the 
northeast of W58 (i.e., east of K3-50C). In this conjuncture, one would expect Sh 2-99 and K3- 
50D to be relatively evolved regions and K3-50C to be the least evolved source as we move from 
southwest to northeast. On a large scale, it seems to be consistent with the scenario of sequential 
star formation because the probable age of K3-50D is ~ 1 - 2 Myr, whereas the expected age 
of K3-50C1 is of ~ few x 10^ yr. This is further supported by the presence of more intense 
molecular material observed towards K3-50C than Sh 2-99 (Israel 1980), along with the average 
extinction gradient of 4.2 to 97 mag from K3-50D to K3-50C1. However, it is not known whether 
the different condensations of Sh 2-100 region were in different stages of evolution before the 
passage of ionization front generated by the stellar winds of WR 1 3 1 . It is also difficult to anticipate 
the projection/ orientation effect on the location of the WR star with respect to the Sh 2-100 region. 
The initial mass of the cloud associated with the individual region is also an important parameter 
as it governs evolution of the cloud in presence of external radiation pressure. Numerical models 
predict that if the pressure due to external ionized gas is too small to compress a massive molecular 
gas, the radiation field will not have any dynamic effect on the cloud (Bertoldi & McKee 1990; 
Lefloch & Lazareff 1994). At present the sequential star formation around Sh 2-100 region due to 
the interaction of H i bubble is suggestive. A more complete picture can be revealed with the help 
of high resolution, high sensitivity H i and CO observations covering the entire W58 complex. A 
positional and chronological study of the YSOs associated with the complex (e.g., Lee et al. 2005; 
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Ogura 2006) can support the validity of sequential star formation. Since the region studied here 

is deeply embedded and shows a large variation in extinction, we would like to determine the 
intrinsic properties of the low mass PMS stars identified from the NIR and IRAC CC diagrams 
with follow up deep optical/NIR spectroscopic observations in order to determine their relative 
ages by constructing intrinsic Hertzsprung-Russell diagrams. These observations will allow us to 
have a better understanding of star formation history across the entire region and a much clearer 
picture of the triggered star formation across the W58 complex. 



8. Conclusions 

In this paper we have presented a multi-wavelength study of the stellar contents and physical 
environment of star-forming region (K3-50) in the proximity of diffuse H u region Sh 2-100. K3- 
50 consists of seven H u regions (A, B, CI, C2, D, E and F) of different classes and are thought to 
be a part of W58 molecular cloud complex located at a distance of ~ 8.7 kpc. Our main conclusions 
are as follows: 

1. We have identified the ionizing star of the optical visible nebula K3-50D. Optical spec- 
troscopy suggests that an 04V star is responsible for the ionization. We have estimated a distance 
of ~ 8.6 kpc for the K3-50D region. The approximate distances for the K3-50E and K3-50F regions 
have been estimated from the NIR photometry and radio continuum flux densities following the 
methodology by Comeron & Torra (2001). The distance estimation suggests that both the regions 
are likely the part of the W58 complex. 

2. We identified the probable exciting sources for six H n regions using the / vs. J - H CM 
diagram. The photometric spectral types of the ionizing sources agree well within a subclass with 
that derived from radio observations based on the number of Lyc photons. 

3. We derived physical parameters for all the five H ii regions using 1280 MHz radio contin- 
uum observation. The H ii regions cover a wide range in size (~ 0.2 to 1.6 pc), optical depth (~ 
0.18 to 1.92), rms electron density (~ 740 to 6390 cm"^) and rate of Lyc photons (~ 5.5 x 10"*' s'^ 
to ~ 1.1 X 10"^^ s~^). These physical parameters reflect different evolutionary stages of H n regions. 

4. We estimated the total mass of the cloud associated with one of the H n regions (K3-50A) 
which was found to be ~ 7800 M©. At the center of the cloud core, we detected a possible massive 
YSO with Mi, ~ 42 + 4.4 M©, accreting with an effective envelope accretion rate of Log(M£„v) ~ 
-3.8+ O.2M0yr-i. 

5. Combining IRAC photometry with ground-based NIR observations within an area of ~ 8' 
X 8', we found a total of 150 objects with circumstellar disks (i.e., YSOs), suggesting a site of 
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active star formation. The distribution of these YSOs (Class II, Class I and H - K > 1.8 sources) 
in Sh 2-100 SFR correlates well with the association of gas and dust within the region. 

6. The distribution of YSOs at the periphery of H i shell on a large scale supports the specula- 
tion by Israel (1980) that the star formation activity around K3-50 region might have been induced 
by an expanding spherical bubble which was created by strong stellar winds from older generation 
OB stars of the complex in the south-eastern direction. 

The Sh 2-100 SFR region represents a broad sample of different stages of massive star for- 
mation all in the same cloud, from the cold, dense core at source CI to the diffuse H n region at 
source F. The conclusions drawn in this paper on this cloud should be relevant and applicable to 
further deep study of the region with high resolution and high sensitivity instruments. 
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Appendix A: Selected interesting region 

IRAS-B is situated at ~ 3'. 3 northwest of K3-50A. An analysis of the IRAS fluxes reveals that this 
source has color characteristics (e.g., log (F25 )/log (F12 ) = 2.02 and log (F60 )/log (F12) = 5.06) 
consistent with being an UCH n region (Wood & Churchwell 1989). A prominent A^s-band nebulosity 
seen in the direction of IRAS-B is offset by ~ 15" in the western direction from the position of 
IRAS-B (see Fig. 24). Our NIR images reveal that the ^s-band nebulosity contains few point 
sources and the position of the MSX point source (G070. 3 164+01 .6493) falls on the ^s-band nebulosity. 
The offset of IRAS-B may be due to the large IRAS beam. Therefore, the more accurate posi- 
tion of the IR source associated with IRAS-B can be inferred from the nearby MSX point source 
(G070.3 164+01. 6493), since the MSX survey has global absolute astrometric accuracy of 1".9. 
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The MSX CC plot (see Fig. 18) for the source suggests IRAS-B to be a massive YSO. Since there 
is no study available in the literature for IRAS-B, the distance to the source is unknown, which 
puts a limit to support in favor of a physical connection of this source with Sh 2-100, rather than a 
chance of alignment along the line of sight. However, if we assume a distance 8.7 kpc, our analy- 
sis leads to many useful conclusions. We acknowledge that this assumption is highly speculative, 
hence the results should be taken with caution. Analysis of NIR CM and CC diagrams suggest that 
the most luminous member associated with IRAS-B is a ZAMS star of spectral type ~ 07-08 with 
NIR excess. Since the source shows NIR excess, the spectral class assigned on the basis of CM 
diagram could have been estimated as an early type rather than its real spectral type, but this still 
disagrees with the ZAMS spectral type (earlier than 03) derived from IRAS point source fluxes 
using the equation of Casoli et al. (1986). This discrepancy could be due to large IRAS beam 
and/or presence of low mass stars, which can contribute significantly to the radiation emitted in 
the FIR. If we assume that an 08 star is present in the region, the absence of radio emission (up to 
~ 5 mJy) at 1280 MHz indicates that (i) either the ionizing star is heavily enshrouded by the dust 
so that all its ionizing radiation is absorbed or (ii) the H u region is suppressed by the accretion of 
matter (see, e.g., Yorke 1986) or (iii) the source must be a ZAMS star of spectral type later than BO. 
The spectral type BO has been estimated from the Lyc photons derived using ~ 5 mJy flux (Bcr flux 
of the map at 1280 MHz) and the calibration table of Panagia (1973), for an electron temperature 
of 10"^ K. In order to get a qualitative estimation of the evolutionary stage of the source, we have 
fitted the models by Robitaille et al. (2006), which are already described in section 6.2.3. The SED 
is shown in Fig. 24. Though, the number of data points used for the SED fitting are quite less, 
nevertheless as stated by Robitaille et al (2007), some of the parameters can be constrained more 
narrowly than the others, depending upon the available fluxes. The parameters computed from 
the models as described in section 6.2.3, imply a source of mass M^, = 10.7 +1.8 M© and total 
luminosity LogL^ = 3.4 + 0.2 L© lying behind 6.1 + 0.2 mag of visual extinction. In the absence 
of sub-mm and mm data points, the resulting uncertainties in other parameters may be large, and 
therefore not quoted here. The mass of the object suggests a star of spectral type B2-B3V, but the 
lack of radio emission around the star at such a low visual extinction requires follow up high sen- 
sitivity, high frequency radio observations, along with ^s-band spectroscopy to obtain the accurate 
nature of the massive star and its distance. 
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Table 1 : Log of observational data. 



Date (UT) 


Filter 


No of 


Exposure (sec) 


Telescope 






frames 


per frame 




2006 Oct 24 


U 


3 


1000 


ST 


2006 Oct 25 


U 


3 


300 


ST 


2006 Oct 24 


B 


3 


900 


ST 


2006 Oct 25 


B 


3 


300 


ST 


2006 Oct 24 


V 


3 


600 


ST 


2006 Oct 25 


V 


3 


100 


ST 


2007 Oct 29 


Gr7 


1 


900 


HCT 


2007 June 18 


Gr7 


2 


900 


HCT 


2008 July 21 


J 


120 


10 


IRSF 


2008 July 21 


H 


120 


10 


IRSF 


2008 July 21 


Ks 


120 


10 


IRSF 


2007 June 18 


Ha 


1 


300 


HCT 


2007 June 18 


R 


1 


300 


HCT 


2007 June 18 


[SIX] 


1 


300 


HCT 


2004 Oct 10 


Bry 


21 


90 


HCT 


2004 Oct 10 


^cont 


21 


20 


HCT 



Table 2: Details of the GMRT radio continuum observations. 



Rest freq. (MHz) 


1280 


610 


Dateof Obs. 


25 Jan 2004 


07 Apr 2004 


Phase center 


«2ooo = 20''01'"45.50^ 


0^2000 = 20''01'"45.50^ 




Siom = 33''32'"43.02^ 


^2000 = 33''32'"43.02^ 


Flux calibrator 


3C48 


3C286 


Phase calibrator 


2052+365 


1924+334 


On source int. time 


l(hr) 


3(hr) 


Cont. Bandwidth 


16 (MHz) 


16 (MHz) 


Synth, beam 


3".3 X 3".0 


5".5 X 4".8 


rms noise 


2.2 mJy/beam 


0.9 mJy/beam 
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Table 3: Observed integrated flux and angular size for associated H n regions 



Frequency 






1280 MHz 


1280 MHz 


Source 


Peak Position 


Int. Flux 


Angular Size 




a(J2000) 


5(J2000) 


Jy 




A 


20:01:45.57 


+33:32:42.86 


0.387 


7".l x4".6 


B 


20:01:48.01 


+33:33:03.21 


1.680 


36".4 X 33".3 


CI 


20:01:53.87 


+33:34:15.52 


0.156 


8". 2 X 5". 3 


C2 


20:01:55.08 


+33:34:16.44 


0.093 


5". 8 X 4". 2 


D 


20:01:47.67 


+33:31:50.61 


1.880 


42". 9 X 37". 8 



Table 4: Derived parameters for the H n regions from radio observations 
Source Tc rms EM log Niyc ZAMSj 

(10^ cm"^) (10^ cm-6 pc) (ph. s'^) Spectral type 



A 


1.92 


6.39 


9.84 


48.36 


08 


B 


0.22 


0.88 


1.14 


49.00 


06 


CI 


1.38 


5.04 


7.07 


48.34 


08 


C2 


0.84 


4.89 


4.28 


47.74 


09.5 


D 


0.18 


0.74 


0.95 


49.04 


06.5 



t From Panagia (1973) 



Table 5: Spectral types of the ionizing sources of the H n regions 



Spectral Type 


Source 


NIR estimation t 


Radio estimation (Ref.) 


A 


>05 


05.5 (1) 


B 


05 


05.5 (1) 


CI 




09(1) 


C2 


B0-O9.5 


09.5 (1) 


D 


05 


06(1) 


E 


B0.5 


08 (2) 


F 


BO 


09.5 (2) 



t From this work 

(1) De Free et al. 1994; (2) Israel (1976) 
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Table 6: Observed and dereddened line intensities derived from the optical spectrum of K3-50A 

nebula 



A (A) 


Lines 


Obs. Flux t 


Dered. Flux ^ 


3726, 9 


[On] 


0.299 


3.121 


4340 


Hr 


0.189 


0.591 


4861 




1.000 


1.000 


4959 


[Om] 


0.502 


0.416 


5007 


[Om] 


1.814 


1.378 


5755 


[Nn] 


0.059 


0.016 


5876 


He I 


0.412 


0.094 


6548 


[Nn] 


1.796 


0.208 


6563 


Ha 


24.79 


2.839 


6583 


[Nn] 


6.74 


0.758 


6678 


He I 


0.281 


0.029 


6717 


[Sn] 


0.502 


0.050 


6731 


[Sn] 


0.932 


0.091 


7065 


He I 


0.457 


0.032 


7136 


[Ar hi] 


0.986 


0.065 


7319, 20 


[0 n] 


0.350 


0.019 


7330,31 


[On] 


0.733 


0.040 



t in the unit of = 2.210 X 10" ergs cm'^ s'^ 



$ in the unit of % = 3.203 x 10"" ergs cm'^ s"^ 
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35 ' 



34' 



33 ' 



32 ' 



31 ' 



3Q ' 



+33 29 ' 




2Q-'^02™QQ' 55 



Fig. 1.— POSS II i?-band image of the region around Sh 2-100. The names (A, B, CI, C2, D, 
E and F) represent the different associated H n regions namely K3-50A, K3-50B, K3-50C1, K3- 
50C2, K3-50D, K3-50E and K3-50F, respectively, following the nomenclature by Israel (1976). 
The RA and DEC coordinates are in J2000.0 epoch. 



-42- 




Fig. 2.— Comparison of IRSF/SIRIUS JHK^ photometry with the 2MASS data for the common 
sources. The difference A (IRSF - 2MASS) in mag is plotted as a function of / mag. 



Fig. 3. — JHK^ color-composite image of the Sh 2-100 star-forming region (7, blue; H, green; and 
, red). Locations of the different regions are also marked. The rectangular box shown by white 
Unes to the left represents the region selected as relative control field (see §5.1). The enhanced 
green color seen below the region D is the artifact of the //-band due to a bright star. North is up 
and east is to the left. 
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Fig. 4. — IRSF ^s-band image in logarithmic scale overlaid by the GMRT radio continuum con- 
tours at 610 MHz. The contour levels are at 0.88 x (3, 5, 8, 12, 17, 22, 28, 35, 43, 51, 70, 81, 
93, 102) mJy/beam, where ~ 0.88 mJy/beam is the rms noise in the map at the resolution of ~ 
5". 5 X 4". 8. Individual H ii regions are also marked in the figure. The labelled axes are in J2000 
coordinates. 



-45- 




20 01 5E 54 52 50 48 4S 44 

RIGHT ASCENSION (J2CX»} 



Fig. 5. — High resolution radio continuum map at 1280 MHz for the region around Sh 2-100 with 
a resolution of 3".3 x 3".0. The contour levels are at 2.2 x (3, 4, 6, 8, 12, 16, 20, 24, 30, 36) 
mJy/beam, where ~ 2.2 mJy/beam is the rms noise in the map. The labelled axes are in J2000 
coordinates. Individual H n regions are also marked. 




Fig. 6. — (a) NIR CM diagram of the Sh 2-100 star- forming region, showing the ZAMS locus (thin 
vertical solid line) reddened by Ay = 0, 10, 20, and 30 mag, respectively, and supergiant locus 
(thick dashed line) at zero reddening. The thick solid lines represent the loci of giants reddened 
by Ay = and Ay = 3.3 mag, respectively. All the loci are drawn at the distance of 8.7 kpc. 
Slanting horizontal lines represent the standard reddening vector with a length of Ay = 30 mag, 
drawn from the MS locus corresponding to different spectral types, (b) NIR CM diagram for the 
relative control region, (c) NIR CC diagram of the Sh 2-100 region for the sources detected in the 
JHKs bands. The continuous line marks the locus of MS and the thick dashed line is the locus 
of giant stars (see inset box). The three parallel dashed lines represent the reddening vectors with 
crosses representing a visual extinction of Ay = 5 mag. The locus of CTTS is also shown by a 
dotted line. The horizontal dashed-dotted line represents the J - H color cut-off to distinguish 
NIR-excess stars from probable other kind of sources (see the text), (d) NIR CC diagram for the 
relative control region. 
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Fig. 7. — (a) & (b). IRAC CC diagrams to identify star-forming galaxies, using the method 
described in Gutermuth et al. (2008). The zones where these extragalactic sources are located 
in the CC diagrams are marked with solid lines, (c). IRAC CC diagram to identify YSOs. The 
different zones indicate the locations of Class II, Class I and Class sources from Allen et al. 
(2004). (d). The H - K/K - [4.5] CC diagram. The curved line is the MS locus of late M type 
dwarfs (Patten et al. 2006). The long-dashed line corresponds to the reddening vector (from 
Flaherty et al. 2007). The dashed-dotted horizontal line represents the H - K color cut off for the 
selection of YSOs (see the text). 
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H-K J-H 

Fig. 8. — (a) CC diagram for the sources within 0.5 pc radius of the H n regions. The giant, dwarf, 
CTT loci and the reddening vectors are plotted in a similar manner as in Fig. 6c. (h) J - H vs. / 
CM diagram for the same sources with their Q values. The nearly vertical solid lines represent the 
ZAMS loci at a distance of 8.7 kpc reddened by Ay = 0, 10 and 20 mag. The slanting lines show 
the reddening vectors corresponding to ZAMS spectral type. 
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Fig. 9. — Spatial distribution of the sources found within 0.5 pc radius of each H ii region su- 
perposed on the ^s-band image. The sources are labeled with numbers corresponding to the 
H II regions. North is up and east is to the left. 





Fig. 10. — Images of K3-50D and K3-50A (clockwise from the top left panel) : (a) The continuum- 
subtracted image of Hor (6563 A) + [N ii] (6548 A and 6584 A), (b) Bias and flat-field corrected 
image of [S n] at 6717A+ 6731 A (c) IRSF ^g-band image with two bright sources (marked with 
crosses) inside the semi-circular shell (d) Continuum-subtracted Bry image, marked with a circle 
of ~ 50" radius around central bright star. 
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Fig. 11. — (a) The blue part of the spectrum of the most luminous star (V ~ 13.31 mag) of K3-50D, 
displaying characteristics similar to ~ 04 MS member, (b) The (U - B)/(B - V) CC diagram for 
the stars lying within 50" of K3-50D region. The continuous curve represents intrinsic ZAMS by 
Schmidt-Kaler (1982), shifted along the reddening vector of 0.72 slope for E(B - V)mm = 1-36 
mag. 




Fig. 12. — The observed flux-calibrated optical spectrum of the compact nebula associated with 
K3-50A. 
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Fig. 13. — Contour map of spatial distribution of cold dust emission at 850 //m around K3-50A 
region. The contour levels are drawn above 3cr, where cr is the rms noise (~ 2 mJy) in the map 
and the FWHM of the symmetric 2-D Gaussian beam is 14".0. The labelled axes are in J2000 
coordinates. 
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Fig. 14. — SED for the luminous source (Al) at the center of K3-50A region. The filled circles 
represent flux values while the triangles denote upper limits of the flux values. The solid lines 
show all the models that fit the data reasonably with;^^ - Xtest^ < 3. The dashed line shows the 
SED of the stellar photosphere in the best-fitting model. 
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Fig. 15. — The blue part of the optical spectrum of the star B3 of K3-50B region. 
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Fig. 16. — Low resolution (~ 44" x 27") radio continuum map at 1280 MHz for the regions K3- 
50E and K3-50F. The contour levels are at 5.0 x (3, 4, 5, 6, 7, 9, 13, 16, 20) mJy/beam, where ~ 
5.0 mJy/beam is the rms noise in the map. The labelled axes are in J2000 coordinates. 
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Fig. 17. — The spatial distributions of Class II (asterisks), Class I (filled triangles), sources detected 
only in H and bands with H-K>1.% (filled circles), YSOs detected from H - K/K - [4.5] CC 
diagram (open circles) and MSX point sources (crosses), overlaid on ^s-band image (for color plot 
see online electronic version). The thin contours represent ^^CO emission from Israel (1980) at a 
resolution of 2'. 3 integrated over a velocity range of 13 km s~^ The thin white contours represent 
8 /im emission from Spitzer. The contours start at 54 MJy sr"' (~ 3cr, where cr is the rms noise of 
the map) and the innermost contour corresponds to a value of 450 MJy sr"^. North is up and east 
is to the left. 
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Fig. 18. — CC plot for the MSX point sources around Sh 2-100 region. The limits derived by 
Lumsden et al. (2002) for massive YSOs and H ii regions (dashed lines) are mentioned by their 
areas. The numbers represent the MSX point sources associated with different H ii regions and are 
mentioned in the upper left comer of the figure. 




Fig. 19.— 



SEDS of bright NIR sources within individual H n regions having IRAC counterparts. 
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Fig. 20. — The intrinsic V/B - V CM diagram for the ionizing stars of K3-50D, K3-50E and K3- 
50F regions. The solid line shows the theoretical isochrones taken from Girardi et al. (2002) for 1 
Myr at a distance of 8.7 kpc. 
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Fig. 21. — CM diagram for YSO candidates in Sh 2-100 SFR. The symbols are same as in Fig. 
17. The PMS isochrone of 0.5 Myr (solid curve) and 1 Myr (dotted curve) from Siess et al. (2000) 
are drawn at a distance of 8.7 kpc and zero reddening. The reddening vectors corresponding to 
0.3Mo and AM o, are also shown in the figure. 
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Fig. 22. — 2MASS ^s-band image overlaid with the warm dust contours at 25 yum (black in color) 
and cold dust contours at 60 jum (white in color) taken from HIRES processed images of IRAS. 
The 25 yum contour levels are at 1.5, 2, 3, 4, 5, 10, 20, 40, 60, 80, 90 % of the peak value 43208 
MJy sr"' and the 60 //m contour levels are at 1.5, 2, 3, 5, 7, 10, 15, 20, 25, 40, 60, 80, 90% of the 
peak value 44942 MJy sr ^ Locations of the individual H ii regions are also marked. The labelled 
axes are in J2000 coordinates. 
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Fig. 23. — MSX A-band image marked with area of H i emission (thick solid line) taken from 
Israel (1980). The crosses represent the positions of different sources in the W58 cloud complex 
(see the text). North is up and east is to the left. 
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Fig. 24. — Left: JHK color-composite image around IRAS-B region. Right: SED of the luminous 
source associated with IRAS-B. The filled circles represent NIR and Spitzer flux values, whereas 
the filled triangles denote upper limit of the MSX flux values. The solid lines show all models 
that fit the data reasonably with - Xtest^ < 3 . The dashed line shows the SED of the stellar 
photosphere in the best-fitting model. 



